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Abstract. Three configurations of two perpendicular disks in M 3 are 
examined, the first in which the disks share centers and the other two in which 
the disks touch at precisely one point. Volume, surface area and mean width 
calculations dominate the discussion. Integrated mean curvature also appears 
as an indirect way to compute mean width. 

Our investigation begins with a theoretical question about experimental data. 
Example 1 is the convex hull of the following two orthogonal disks in R 3 : 

{(x,y,z) : x 2 + y 2 < 1 & z = 0} and {(x, y, z) : x 2 + z 2 < 1 & y = 0} . 

We can numerically evaluate the volume VL, surface area AR and mean width MW 
of the corresponding solid domain in Figure 1 using [lj: 

VLt ps 2.666, ARi ps 10.28, MW X ps 1.869. 

Example 2 is the convex hull of the two disks: 

{(x, y, z) : x 2 + y 2 < 1 & z = — 1} and {(x, y, z) : x 2 + z 2 < 1 & y = 1} 

with corresponding solid domain in Figure 2 and 

VL 2 ~ 3.141, AR 2 ps 13.92, MW 2 ps 2.277. 

Example 3 is the convex hull of the two disks: 

{{x,y,z) : x 2 + (y + l) 2 < 1 & z = 0} and {(x, y, z) : (y - l) 2 + < 1 & x = 0} 

with corresponding solid domain in Figure 3 and 

VL 3 ps 3.627, AR 3 ps 15.97, MVT 3 ps 2.645. 

Of the nine constants, just two (VLi = 8/3 and V L 2 = vr) are readily identifiable. 
What are exact closed-form expressions for the remaining constants? 
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Given Q to be a convex body in M 3 , a width is the distance between a pair of 
parallel f2-supporting planes. Every unit vector v G M 3 determines a unique such pair 
of planes orthogonal to v and hence a width w(v). Let v be uniformly distributed 
on the unit sphere S 2 C 1R 3 . Then w is a random variable and its average value is 
the mean width of Q. Three numerical characteristics of Q - volume, surface area 
and mean width - are central to our study. These quantities, along with the Euler 
characteristic, form a basis of the space of all additive continuous measures that are 
invariant under rigid motions in R 3 . 

"The mean width is a new measure on three-dimensional solids that enjoys equal 
rights with volume and surface area" j2], hence much of this paper is devoted to 
computing MW for our three examples. What we call the direct approach is based 
on the definition of MW; what we call the indirect approach utilizes a connection 
between MW and integrated mean curvature (often called "integral" or "total" mean 
curvature). This connection is suggested in the materials science [3J H] and astro- 
physics literature [5J E]; the closest claim to a proof appears in [7], based chiefly on 
[8]. Our paper therefore also serves to confirm the validity of the indirect approach 
for certain non-polyhedral test casesQ 

1. Example 1 

The boundary dQ of the convex hull Q here is trivially given by the surface 



z = ± (VT 

over the planar region x 2 + y 2 < 1. Let 



x 2 — \y\ 



(p(x, y) = VI — x 2 — y, x 2 + y 2 <lky>0 
then (p x , ip y , (p xx , <f , ip yy denote first /second-order partial derivatives of cp and 



1 vT^x 



2 



VLx = aJ J <p(x,y) dydx = ~, 



-1 o 



1 Vl-x 2 



AR X = JdS = 4:J J yjl + ipl + ip 2 y dy dx = 2(2 + tt) 

an -l o 

= 10.2831853071795864769252867.... 



page 513 of [7], mean curvature K. is defined as the average of the two principal curvatures, 
but this is inconsistent with [8], which takes JC to be the sum. We follow [8], defining 2H = K,. Our 
formula correctly gives MW = (£ + irr)/2 for a right circular cylinder of length I, radius r [9l HOj. 
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1.1. Indirect Approach. Let dQ + denote the upper portion of dfl and dQ~ 

denote the lower portion. On the one hand, the mean curvature of dQ + is 



H(x,y) = 



(! + Vl) Vyy ~ ^VxVyVxy + (l + ¥>J) ¥z 



2(1 + ^ + ^) 3/2 (2-^)3/2 
over the open region x 2 + y 2 < 1 & y > 0. It follows that 



H dS = H(x, y) + + dy dx = - _ ^ ^—^ dy dx = {2 _J ) ^ l _ x2 - 

On the other hand, the exterior dihedral angle on the semicircular edge x 2 + y 2 = 1 
& y > is 

1 



a = 2 arccos 



,v / 2^" 2 

because the unit exterior normal vector to dfl + is 



and the unit exterior normal vector to the cylinder x 2 + y 2 = 1 is (x, y,0). The 
dot product of the two vectors is l/y/2 — x 2 ; we multiply the angle by two since 
the dihedral angle between dQ + and dQ~ is twice the preceding angle. In terms of 
arclength s = 9, < 9 < ir, we have 



ads = 2 arccos = d9. 

\V2-cos 2 9J 

The surface <9f2 is piecewise continuously different iable and has n = 4 smooth edges 
Ej with (non-constant) dihedral angles a,-, 1 < j < n. From the general formula 



MW = ±-J H dS + l-<rJa J ds, 



we deduce that 



1 y/T— X 2 7T 

MWi = — / ; + — / arccos - d9 

2nJ J (2-x 2 )VT Zr ^ 2 W VV2-cos 2 0/ 



-1 o 



-- ^v^ln (-1 + v 7 ^ -4Li 2 (-1 + v 7 ^) +4Li 2 (l - v 7 ^) 



= 1.8697727582861870379136441... 
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where Li 2 is the dilogarithm 

X '\n(l-t) 



i 1 



t 

k=l 



dt, \x\ < 1. 



1.2. Direct Approach. Consider the portion of dQ in the first octant only. In 
this octant, an ^-supporting plane P d : 

—x H — -y H — -z = 1 (with coefficients a > 0,b > 0, c > and scaling factor d > 0) 
a a a 



has an associated line Lt„ 



a b 

-x + -y = lkz = 
d d 

in the xy-plane and an associated line L d xz : 

-x + -z = 1 k y = 
a a 

in the xz-plane. Assume WLOG that a 2 + b 2 + c 2 = 1. The distance of P d from the 
origin O is d Also, 

distance of L d , from O is — = and distance of Li, from O is 



The largest d such that 



is thus 



supports x 2 + y 2 = 1 or L^ 2 supports x 2 + z 2 = 1 



d = max | V a 2 + b 2 , V a 2 + c 2 j . 



Let 

a = cos sin 0, 6 = sin 9 sin 0, c = cos 

where < < ir/2, < < ir/2. To ensure uniformity, think of (9, 0) as possessing 
joint density - sin 0. We have 



7r/2 7r/2 



MW 7 ! = 2 J J max |sin 0, -y/ cos 2 sin 2 + cos 2 1 — sin <i0 
o o 

tt/2 tt/2 

sin 2 <pd<pd9 = 1.8697727582861870379136441... 

7T 
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where 

sin 9 

£(0) = arccos 



V^-cos 2 ^ 

is the required solution (for in terms of 6) of the equation 



sin = \J cos 2 9 sin 2 + cos 2 0. 



2. Example 2 

The curved portions of the boundary dVt of the convex hull Q here are given by 

z = (p(x, y) and z = i/)(x, y), where 

(p(x, y) = y/l — x 2 + y — 1, |x| < 1, \y\ < 1 & y > — VI — a; 2 , 

4>(x,y) = —y/1 — x 2 + y — I, \x\ < 1, |y| < 1 & y > VI — a; 2 - 

The flat portions of dQ are the two disks, one of which is given by z = —1 over 
x 2 + y 2 < 1. These facts contribute to the following: 

l v 7 ! 71 ^ i i 

VL 2 = J J (<p{x,y) + 1) dydx + J J (ip(x, y) - ip(x, y)) dy dx = ir, 

ii ii 

AR 2 = 2tt + J J JT+tfTrfdydx + J J ^l + ^ 2 x + ^ 2 y dydx 
= 13.9235808852350105127348109... 



where 



tt/2 1 

E(n) = J v 7 ! -I2sm(6) 2 d6 = J < 



l -^dt 



1-t 2 





is the complete elliptic integral of the second kind. 



2.1. Indirect Approach. Let <9f2 + denote the curved portion of <9f2 prescribed 
by if and dQ~ denote the curved portion prescribed by ip. We have 

( 1 + <Pl)<Pyy- 2 <Px<Py<Pxy+ { l + ^y) Vxx _ 1 



2(l + ^ + ^) 3/2 ~ (2-x 2 ) 3 / 2 

(1 + V') ^yy ~ ^x^y^xy + (l + tf) ^ a 



2(1+^ + ^) 



2\3/2 
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everywhere and hence 

dy dx 



HdS 



(2 - x 2 ) s/i^x 2 
as previously. It follows that 

fHdS=[ j dydQ ; + / / dyd * =V2«. 

J J J (2-x 2 )VT^ J J {2-x 2 )Vl^x 2 

an -i_yi^2 -171^2 

Let e denote the circular edge x 2 + y 2 = lhz = — 1. Clearly dQ + C\e is the semicircle 
with y = —y/l — x 2 whereas d£l~ He is the semicircle with y = y/l — x 2 . The exterior 
dihedral angle on e is 

( v 

a = arccos . 
because the unit exterior normal vector to dQ + , dQ~ is 



' (- Vx , -<p y , 1) = J LJL (- 7 J==, -1,1) = - 7 =L= (x, y, -y) , 



X 2 



(*..*.. - 1 ) = vS (71%' _1 ) = Trb? ^ _!/) 



1 /, , _n \\ — X 2 ( X 

respectively and the unit exterior normal vector to the horizontal disk is (0,0, — 1). 
The dot product of the two vectors is y/y/2 — x 2 . An identical argument applies for 
the circular edge x 2 + z 2 = 1 & y = 1. In terms of arclength s = 9, < 9 < 2n, we 
obtain 



2tt 



ljads = 2j 



e 

which leads to the conclusion that 



arccos I = I dv = 2n 

V2-COS 2 0, 



MW 2 = — (V2tt\ + (2tt 2 ) = - (V2 + n) = 2.2779031079814441436321660.... 
2n V / An 2 V / 

2.2. Direct Approach. Consider the curved portion of dQ in the halfspace x > 
only. In this halfspace, an ^-supporting plane P d : 

-x + -v + -z = 1 (with coefficients a > 0, b, c and scaling factor d > 0) 
d d d 

has associated lines L^ y , : 

a 6 d + c a c d — b B 

-x + -y = — — & z = -1, -x + -2 = — — & y = 1. 

a a a a a a 
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Assume WLOG that a 2 + b 2 + c 2 = 1. The distance of P d from the origin O is d. 
Also, 

distance of Lt , from O z is = - and distance of L d from 0„ is 



where O z = (0,0, —1) and O y = (0, 1,0). The largest d such that one of the unit 
circles is supported is thus 



d = max | V a 2 + b 2 — c, V a 2 + c 2 + 6 j . 



We introduce spherical coordinates as before, but with — 7r/2 < < tt/2, < <p < tt 
instead. To ensure uniformity, think of (9, (ft) as possessing joint density 7^rsin0. 
We have 

7r/2 7T 

MW2 = 2 I I max I sin <p — cos <p, \J cos 2 sin 2 <p + cos 2 + sin sin 1 — sin <p d<p dQ 



-tt/2 

sin d> — cos d>) smd>dd>d6 = 2.2779031079814441436321660... 



2 

7T 

where 

£(0) = vr/2 + arctan(sin0) 
is the required solution (for in terms of 0) of the equation 



sin cj) — cos (f) = y cos 2 sin + cos 2 + sin sin 0. 

3. Example 2 (Again) 

Vinzant, using techniques in her thesis [11 j . computed that dQ is given implicitly by 
the equation 

= —x + 2y + x y — 3y + y — 2z — x z + 6yz 

+a; 2 ?/z - 5y 2 z + y 3 2; - 3z 2 + 5yz 2 - 2y 2 z 2 - z 3 + y^ 3 
= {y-l){z + l)(x 2 -2y + y 2 + 2z-2yz + z 2 ) 
= (y-l)(z + l)(z-if(x,y))(z-i/;(x,y)) 

verifying what we already know. She additionally gave an elegant parametric repre- 
sentation of the curved portion in x > 0: 



x 



v 7 ! — v 2 , y = 1 — u + uv, z = —u — v + uv, < -u < 1 & — 1 < v < 1 
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which deserves further attention. In the following, we reproduce our results from 
the preceding section. The purpose in doing so is not to torture the reader, but 
rather to set the stage for Example 3 (for which a parametric representation is the 
only workable method available.) The Jacobian determinant^! 



d{x,y) 
d{u, v) 



Xu X v 



(-1 + v)v 



allows us to evaluate 



Defining 



we have 



Defining 



i i 



VL 2 = 2 / / (( — u -v + uv) + l) ^ 1+t '^ dv du = 7T. 

V 1 — f 2 



o -1 



E 



F \Xui Uui Zu) ' y^V) Vvi %v) 



1 1 



AR 2 



dS = 2n + 2 / VEG - F 2 dv du 



mi 



2tt + 2 



o -1 



i i 



(1 -v)(l + v 2 ) 



1 + v 



dv du 



o -1 



2^ + 2^(1)). 



\{x u , Vu, z u ) x [x V j y Vj z v ) | 



1 -v 2 



M = -- ((x u , y u , z u ) ■ Af v + (x v , y v , z v ) ■ Af u ) 



2 The fact that this has indefinite sign doesn't affect the volume calculation. 
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we have 



EN-2FM + GL 
HdS = 1 2{EG-F 2 ) dS 



on an 

1 i 



o -1 

i i 



EN - 2FM + GL 
2 {EG - F 2 ) 



VEG - F 2 dv du 



1 ./a-* 1 



o -1 

i i 



(l+v 2 ) 3/2 V l + v 



= 2 / / -J- -dvdu = V2n. 

J J 1 + v 2 V 1 + v 

-1 

The semicircular edge x 2 + z 2 = l,y = l&zx>0 corresponds to m = & — 1 < v < 1. 
Call this e. The exterior dihedral angle is 

a = arccos (J\f ■ (0, 1, 0)) = arccos 

and arclength s satisfies 



Vl + v 



Consequently 



and 



ds = \J x\ + z 2 dv = —j^=^= dv. 

Vl — v 2 



i 

4 / a ds = 4 / , arccos ( , ) dv = 2n 2 

VT^v 2 KVTTT 



as was to be shown. 

4. Example 3 

Vinzant, using techniques in her thesis [TT] . computed that dQ here is given implicitly 
by the equation 



= -4x 4 + 8x 6 - 16x 2 y + 36x 4 y - 16y 2 + A0x 2 y 2 + 15x 4 y 2 + 36x 2 y 3 + 8y A 
+6xV - y 6 - 8x 2 z 2 + 24xV + IQyz 2 + A0y 2 z 2 - 18x 2 y 2 z 2 - 36y 3 z 2 
+6y 4 z 2 - 4z 4 + 24xV - 36yz A + lhy 2 z A + 8z 6 . 
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One could solve for this cubic (in z 2 ) and proceed as earlier, laboring against the 
weight of complicated expressions. We prefer, however, to exploit another of her 
elegant parametric representations: 

x = U j-v{2 + v), y= — , z = _ 

for < u < 1 & —2 < v < —2/3. The Jacobian determinant 



d(x, y) I —v 2 + v + 6uv + Quv 2 



d(u,v) \2 + v (l + 2v) 2 

allows us to evaluate 

1 -2/3 

A f f (-l + u)Jv(2 + 3v) I -v 2 + v + Quv + Quv 2 , , 

KL 3 = 4/ / — W- - — - dvdu 

J J l + 2v \ 2 + v (l + 2v) 2 

-2 

= -J = 3.6275987284684357011881565.... 

Likewise, 

AR 3 = A I I VEG -F 2 dvdu = 4 / / "-rcTu^-rucc , .L-r^-ruc , 

J J J J (l + 2v) 2 V (2 + t;)(2 + 3*;) 

0-2 0-2 

-2/3 



1 —2/3 1 —2/3 

2 + v + 6m; + 6uv 2 / 1 + 2v + 3v 2 



2 + 4w + 3v 2 / l + 2w + 3v 2 , 

—, w x rfw = 15.9716335277272626753537899.. 

(l + 2i;) 2 V (2 + i;)(2 + 3v) 

-2 V 

Gosper & Bickford [12J conjectured that Ai? 3 = 4cx/r, where 



-8HE f-^ (l7-56i\/2^ + 9^/-17 + 56r\/2£ ^ (l7 + 56n/2^ 
+ 108v / 27<" - 72^^ (l7 + 56iV2^ 
- (56^ + 32^) n (7 - 4zV2) , ^- (l7 - 56iV2 H 
+ ^56i + 32^) n Q (V - 4^) , ^ (l7 - 56iV2 
+72* n ^ (7 + 4iv^) , ^ (l7 + 56^) J 
-12% n Q (7 + 4iV2) , (l7 + 56^) 
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and r = 18 (—i + 2y/2), where i is the imaginary unit, E(fj,) was defined earlier, 

n/2 1 

KU)= I — 1 =d0= [ - 1 rft 
7 Vl-Zisin^) 2 7 v/(l -t 2 )(l - /it 2 ) 

is the complete elliptic integral of the first kind and 

tt/2 1 

n(z/,u)= / 1 ; =d9= [ 1 dt 

J (l-z/sin(#) 2 ) y/l-fj,sm(e) 2 J (l-z/t 2 )v/(l-t 2 )(l-/it 2 ) 

is the complete elliptic integral of the third kind. A simplification of a would be 
good to see someday. 

4.1. Indirect Approach. As before, 



= (x u ,yu,z u ) x (x v ,y v ,z v ) _( J -v(2 + v) 1 + v j v(2 + 3v) 



)\ \ V 1 + 2v + 3v 2 ' y/i + 2v + 3v 2 ' V 1 + 2v + 3v 2 

and 

1 -2/3 

3t;(l + 2t;) 2 



-2 

1 -2/3 



-4 



-2 



[l + 2v + 3v 2 fl 2 {2 + v + 6uv + 6uv 2 ) 



2 + v + 6uv + 6uv 2 / l + 2v + 3v 2 , , 

-, \ w r dv du 

(l + 2v) 2 y (2 + v)(2 + 3v) 

-2/3 

= -4 / ; dv 

J 1 + 2v + 3v 2 ^T(2T^)(2T3^) 

= 2^. 

Let £ denote the arc of the semicircle x 2 + (y + I) 2 = 1, z — 0Szx>0 that runs 
counterclockwise from points (0, —2, 0) to (2\/2/3, —2/3, 0); this corresponds to u — 1 
& —2 < v < —2/3. The exterior dihedral angle is 



a = 2 arccos (AT • ^-v{2 + v), 1 + v, o)) = z arccos , 



2v + 3v 2 
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because the unit exterior normal vector to the cylinder x 2 + (y + 1) 2 = 1 is (x, y + 1, 0); 
the arclength s satisfies 



ds= \fxl + y 2 v dv = dv. 

y/-v{2 + V) 

Consequently 

-2/3 

4 ads = 8 / — = arccos ( . ] dv 

J J J-v(2 + v) \Vl + 2v + 3v 2 , 

e -2 v 

= 47rarcsec(3) 

and therefore 

MW 3 = -^lHdS+-lads 



2lX J 7T 

dfl e 

= V2 + arcsec(3) 

= 2.6451729797138697309366179.... 

All nine constants exhibited (at the beginning) possess closed-form expressions, al- 
though the result for AR 3 is partly conjectural. We had expected that there might be 
required "more time to develop the languages, functions, symmetries, etc., to express 
the constants more naturally" [12], but this belief turned out to be overly cautious. 

4.2. Direct Approach. Consider the portion of dQ in the first octant only. In 
this octant, an f2-supporting plane P d : 

a b c 

— x + —y + — z = 1 (with coefficients a > 0, b > 0, c > and scaling factor d > 0) 
a d d 

has associated lines lA,„ L d : 

(lb be 

—x + —y = 1 & z = 0, —y + —z = 1 & x = 0. 

a d d d 

Assume WLOG that a 2 + b 2 + c 2 = 1. The distance of P d from the origin O is d. 
Also, 

distance of Lt„, from 0_ is — ; and distance of L d from 0+ is 



*» " ^2^52 y* + VV r Td 2 ~ 
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where 0_ = (0, —1,0) and + = (0, 1,0). The largest d such that one of the unit 
circles is supported is thus 

d = max | V a 2 + b 2 — b, Vb 2 + c 2 + 6 j . 

We introduce spherical coordinates with O<#<7r/2,O<0<7r/2, obtaining 

7r/2 tt/2 

MW3 — 2 / / max < (1 — sin 9) sin 0, sin 9 sin + -y/ sin 2 # sin 2 + cos 2 i — sin <i0 d# 








sin sin + y sin # sin + cos 2 sin c/0 d9 





ft tt/2 



-J J (1 -sin#)sin 2 0d0d# 



4 





tt/2 tt/2 



+ — y J ^sin sin + y sin 2 9 sin 2 + cos 2 ) smo<7o<76> 

K 

2.6451729797138697309366179... 



where k = arccsc(3) and 



«») = axccosh/i^±|^|], 0<0< K 
2 — 4sm# + 3 sin 9 ' 



is the required solution (for in terms of 9) of the equation 



(1 — sin 9) sin = sin # sin + y sin # sin + cos : 

5. Related Topics 
Dirnbock & Stachel HU studied the convex hull of the two disks: 



Ux,y,z) :x 2 + (y + f) 2 < 1 & ^ = 0} and y, 2) : (y - f f + z 2 < 1 & x = o} 



when 5=1 and Ira [15] studied the same when 5 = \/2. These are intermediate 
cases relative to our Example 3 (for which 5 = 2) and what is essentially Example 1 



Convex Hull of Two Orthogonal Disks 



14 



(for which 5 = 0). The former case, called an oloid, has volume 

V L = If- < 2 + mB9 » 2 W 

3 J (1 + cos 9) Vl + 2cos0 

- I(-^fr9^fri)) 

= 3.0524184684243748566972005... 

and surface area AR = 4n, where E & F are incomplete elliptic integral of the second 
& first kinds respectively: 

<j> sin 4> j 

£(<M = Jy/l-pBm(0)*d9= J J^z^dt, 





sin cfr 



^ I Jl-usm(6) 2 M I 



(of course, F(n/2, •) = if (•)). A feature of the oloid is that each of the two inter- 
locking orthogonal circles C\ c f t & C r i g ht intersects the center of the other. Define 

constants to = arcsin ( v a/2 — 1 ) and 



1 f (2 + v^oostf) 2 (l + cos9 + cos 9) 

sy — / — — dQ 

2V2J (i + v^cos^) 2 + cos 9 + cos 2 9 
= — + E(u, -1) + V2 (U(-V2 - 1, to, -1) + n(v / 2 - 1, cu, -1) - F(u, -1) 
where II is the incomplete elliptic integral of the third kind: 

<j> sin<j> 

n(z/, 0, a) = / d6= dt 

J (l-z/sin(#) 2 ) y/l-(j,sm{6) 2 J (1 - vt 2 ) y/(l - t 2 )(l - fit 2 ) 

The latter case, called a two-circle roller, has volume 

VL = -^=7 = 3.2818194874496894190321933... 
3v 2 
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and surface area 

AR = 87 = 13.9235808852350105127348109.... 

As the name "roller" suggests, the authors of [TH [T2] were raising issues of a physi- 
cal/mechanical nature. We suspect that an algebraic approach based on [TT] could 
supplant some of their technical arguments leading to VL & AR. Exact expressions 
for MW remain open, as far as is known. 

A picture in [16] depicts the circles CWt & C r i g ht but not within our context of 
convex hulls (taut rubber sheets spanning wire frames). The context instead is about 
minimal surfaces (elastic soap films spanning the same), for which mean curvature is 
zero everywhere. Many questions suggest themselves. 

6. Acknowledgements 
Wouter Meeussen's package ConvexHull3D.m was helpful to me in preparing this 
paper [I]. He kindly extended the software functionality at my request. Cynthia 
Vinzant generously computed the parametric representations for both Examples 2 
and 3 based on |TT] . I am also thankful to R. William Gosper, Neil Bickford, Roland 
Roth, Rolf Schneider, Hiroshi Ira, Thinh Le, Qiang Du, Frank Sottile and Tina Mai 
for their helpful correspondence. 

7. Addendum 

Axel Vogt improved upon the conjecture in [12] and successfully obtained 



AR 3 = - 
d 3 



via a change of variables and Maple. It is possible to reduce the complicated expres- 
sion 4a /t further: 

a = -2431E (l7 - 56i\/2) j + 27\/-17 + 56iV2E (l7 + 56iV2^ 
+ 162^^ Qj - 162^^ (l7 + 56iV2^ 
+6 (8i - 7V?) n (J^ (7 - Uy/tj , (l7 - 56iV2j\ 

+2 (lOAi + 71 v^) n Q ( 7 ~ ' ^1 ( 17 ~ 56 ^ 

(-i + 2V2 



t = 54 

using numerics and Mathematica, but Bill Gosper's question (on symbolic transfor- 
mations between elliptic integrals to link these) remains unanswered. 
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Figure 1: Centers coincide 
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Figure 2: Centers on diagonal 
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Figure 3: Centers on bisecting line 
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